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Abstract 

We present a classification scheme for pliase transitions in finite systems like 

atomic and molecular clusters based on the Lee- Yang zeros in the complex 

temperature plane. In the limit of infinite particle numbers the scheme reduces 

to the Ehrenfest definition of phase transitions and gives the right critical 

indices. We apply this classification scheme to Bose-Einstein condensates in 

a harmonic trap as an example of a higher order phase transitions in a finite 

system and to small Ar clusters. 
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Small systems do not exhibit phase transitions. Following Ehrenfest's definition this 
statement is true for almost all small systems. Instead of exhibiting a sharp peak or a 
discontinuity in the specific heat at some well defined critical temperature the specific heat 
shows a more or less smooth hump extending over some finite temperature range. For 
example for the melting of atomic clusters this is commonly interpreted as a temperature 
region where solid and liquid clusters coexist and as a finite-system analogue of a first 
order phase transition. Proykova and Berry |^ interpret a structural transition in TeFg 
clusters as a second order phase transition. A common way to investigate if a transition in 
a finite system is a precursor of a phase transition in the corresponding infinite system is 
to study the particle number dependence of the appropriate thermodynamic potential 0. 
However, this approach will fail for all system types where the nature of the phase transition 
changes with increasing particle number which seems to be the case e.g. for sodium clusters 
or ferrofluid clusters . For this reason a definition of phase transitions for systems with 
a fixed and finite particle number seems to be desirable. The only recommended feature is 
that this definition should reduce to the Ehrenfest definition when applied to infinite systems 
for systems where such limits exist. An approach in this direction has been made by studying 
the topological structure of the n-body phase space and a hypothetical dehnition based on 
the inspection of the shape of the caloric curve |0] . A mathematical more rigid investigation 
giving the sufficient and necessary conditions for the existence of van der Waals-type loops 
has been given by Wales and Doye 0. 

Our ansatz presented in this letter is based on earlier works of Lee and Yang |^ and 
Grossmann et al. who gave a description of phase transitions by analyzing the distri- 



butions of zeros (DOZ's) of the grandcanonical and the canonical partition function 
Z{(3) in the complex temperature plane. For macroscopic systems this analysis merely con- 
tributes a sophisticated view of the thermodynamic behavior of the investigated system. We 
will show that for small systems the DOZ's are able to reveal the thermodynamic secrets 
of small systems in a distinct manner. In the following we restrict our discussion to the 
canonical ensemble and denote complex temperatures hy B = (3 + it where (3 is as usual 



l/k^T [|n . 



In the case of finite systems one must not deal with special considerations 
regarding the thermodynamic limit. We write the canonical partition function 
Z{B) = J dEfl{E) exp{—BE), with the density of states as a product 

ZiB) = ZU>3)Z;^,iB) (1) 

where Zii^{B) describes the limiting behavior of Z{B) for T — oo imposing 
limg^o ^int('B) = 1. In general, Z^^^B) will not depend on the interaction between the 
particles or the particle statistics but it will depend on the external potential imposed. 
E.g. for particles in a rf-dimensional harmonic trap we have Zii^{B) ~ q-^n 
(i-dimensional gas Zi^j^iB) ~ the following we will assume that Zii^{B) has no 

zeros except at B = oo. Then the zeros of Z{B) are the same as those of Zi^tiB). Applying 
the product theorem of Weierstrass [|1^ the canonical partition function can be written as 
a function of the zeros of Z[^t{B) in the complex temperature plane. Because Z{B) is an 
integral function its zeros Bk = B*_^ = Pk + i^k {k G N) are complex conjugated 

Z{8) = Z„„(/J)Z,„(0)exp(eaBlnZ,„.(0)) (2) 
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The free energy F{B) = —^\ia{Z{B)) is analytic, i.e. it has a derivative at every point, 
everywhere in the complex temperature plane except at the zeros of Z{B). If the zeros are 
dense on lines in the complex plane, different phases are represented by different regions of 
holomorphy of F{B) and are separated by these lines in the complex temperature plane. The 
DOZ contains the complete thermodynamic information about the system and all desired 
thermodynamic functions are derivable from it. The calculation of the specific heat Cv{B) 
by standard differentiation yields 

1 1 



Cv{B) = a^{B)-kBB'J2 



_{Bk-BY~^ {Bl-B)\ 
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Zeros of Z{B) are poles of F{B) and Cv{B). As can be seen from eq. the major contri- 
butions to the specific heat come from zeros close to the real axis, and a zero approaching 
the real axis infinitely close causes a divergence in the specific heat. 

In the following we will give a discretized version of the classification scheme of Gross- 
mann et al. To characterize the DOZ close to the real axis let us assume that the zeros 
lie approximately on a straight line. The crossing angle of this line with the imaginary axis 
(see Fig. |lD is then u = tan 7 with 7 = {(32 — Pi) / {t2 — ti) . The crossing point of this line 
with the real axis is given by Pent = Pi — l^i- We define the discrete line density as a 
function of as the average of the inverse distances between Bk and its neighboring zeros 

Guide lined by the fact that the importance of the contribution of a zero to the specific heat 
decreases with increasing r we approximate 0(r) in the region of small r by a simple power 
law 0(r) ~ r". A rough estimate of a considering only the first two zeros yields 

In0(r3) - In0(r2) . . 

a = . (5) 

In Ts — In T2 

Together with ri, the imaginary part of the zero closest to the real axis, the parameters 7 and 
a classify the DOZ. As will be shown below, the parameter ri is the essential parameter to 
classify phase transitions in small systems. For a true phase transition in the Ehrenfest sense 
we have ri 0. For this case it has been shown |10[ that a phase transition is completely 



classified by a and 7. In the case a = and 7 = the specific heat Cv{P) exhibits a 5-peak 
corresponding to a phase transition of first order. For < a < 1 and 7 = (or 7 7^ 0) 
the transition is of second order. A higher order transition occurs for 1 < a and arbitrary 
7. This implies that the classification of phase transitions in finite systems by 7, a, and ti, 
which reflects the finite size effects, is a straightforward extension of the Ehrenfest scheme. 

The imaginary parts Tj of the zeros have a simple straightforward interpretation in the 
quantum mechanical case. By going from real temperatures P = l/{k-BT) to complex tem- 
peratures P + ir/h the quantum mechanical partition function can be written as 



Z{p + ir/h) = Tr (expi-irR/h) exp(-/3H)) (6) 

= (^'can I exp {-trR/h) I ^can) (7) 
= (^can(t = 0) I ^can(t = t)) , 

introducing a "canonical state" \ \l/can) = I ^i) > which is the sum of 

all eigenstates of the system appropriately weighted by the Boltzmann factor. Within this 
picture a zero of the partition function occurs at times Tj where the overlap of a time evoluted 
canonical state and the initial state vanishes. This resembles a correlation time, but some 
care is in order here. The time x; is not connected to a single system, but to an ensemble 
of infinitely many identical systems in a heat bath, with a Boltzmann distribution of initial 
states. Thus, the times Tj are those times after which the whole ensemble loses its memory. 

Equation (^) is nothing but the canonical ensemble average of the time evolution operator 
exp {—irH / K) . Following Boltzmann the ensemble average equals the long time average 



which was proven quantum mechanically by Tasaki [|T^. Therefore Tj indeed resembles 
times for which the long time average of the time evolution operator vanishes. 

The observation of Bose-Einstein condensation in dilute gases of finite number (~ 10'^ — 



10^) of alkali atoms in harmonic traps [M has renewed the interest in this phenomenon 



which has already been predicted by Einstein in 1924. The number of condensed atoms 



in these traps is far away from the thermodynamic limit, raising the interesting question 
how the order of the phase transition changes with an increasing number of atoms in the 
condensate. For this reason we treat the Bose-Einstein condensation in a 3-dimensional 
isotropic harmonic trap {h = u = = m = 1) as an example for the application of the 
classification scheme given above. 

For non-interacting bosons the occupation numbers of an eigenstate \i) and particles 
can be evaluated by a simple recursion [|1^] 

V^{N + 1,B) = -P^^ exp(-i3e.)(r/,(iV, B) + l). (8) 

Since the particle number is a conserved quantity in the canonical ensemble the direct 
calculation of the normalization factor can be omitted by using the relation 



ZN+m ZT=o ew{-Be,){ri,{N, + 1) " ^ ^ 

Since Zpf{B) is an exponentially decreasing function in /? it is a difficult numerical task 
to calculate its zeros directly. Zeros of the partition function are reflected by poles of the 
ground state occupation number 

evaluated at complex temperatures. Fig. 2 displays contour plots of \rjo{N,B)\/N for 40, 
120, and 300 particles. The locations of the zeros of Z{B) (poles of r/o(iV, B)) are indicated 
by the white spots. The separation of the condensed (dark) and the normal (bright) phase is 
conspicuous. The zeros act like " boundary postd^ between both phases. The boundary hne 
between both phases gets more and more pronounced as the number of particles increases 
and the distance between neighboring zeros decreases. Fig. 2 virtually displays how the 
phase transition approaches its thermodynamic limit. We have determined the classification 
parameters for the phase transition by a numerical analysis of the DOZ for up to 400 
particles. The results are given in Fig. 3. The parameter a is constant at about 1.25. 
The small fluctuations are due to numerical errors in the determination of the location of 
the zeros. This value of a indicates a third order phase transition in the 3-dimensional 
harmonic trap. Results for the 2-dimensional systems and other trap geometries, which 
will be published elsewhere in detail, indicate that the order of the phase transition depends 
strongly on the trap geometry. The parameter 7 and the non-integer fraction of a are related 
to the critical indices of the phase transition, e.g. 7 = indicates equal critical indices for 
approaching the critical temperature from the left and from the right. Regarding the finite 
size effects Ti is of major importance. As can be seen in Fig. 3 (b) Ti is approximately 
proportional to 1/iV that the systems of bosons in a 3 dimensional harmonic trap approaches 
a true higher order phase transition linearly with increasing particle number N. 

It appears that the DOZ for Bose-Einstein condensates is rather smooth. As an example 
for a little more complicated situation we calculated the DOZ for small Ar clusters, which 
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have been extensively studied in the past IT^-pO] . Their thermodynamic behavior is governed 



by a hopping process between different isomers and melting pT|-p3|. Many indicators of 
phase transitions in Ar clusters have been investigated, e.g. the specific heat [p4| , p5| , the 
rms bond length fluctuation and the onset of an l//-noise behavior of the potential 
energy in time dependent molecular dynamics simulations . However, for a good reason. 



all these attempts lack a definite classification of the transitions taking place in these clusters. 
Without going into the details of our numerical method which is based on a determination 
of the interaction density of states by extensive Monte Carlo simulations along with an 
optimized data analysis pSf we give here the results for Arg and At^q. Fig. 4 displays 
contour plots of the absolute value of the specific heat cyiB) in the complex temperature 
plane. For Are the poles lie on a straight line at T ~ 15 K and are equally spaced with 
resulting classification parameters a = 0, 7 = 0, and Tih = 0.05 ps. From earlier works 
1 2^ it is well known that at this temperature a hopping transition between the octahedral 



and the bicapped tetrahedral isomer occurs. Our classification scheme now indicates that 
this isomer hopping can be identified as a first order phase transition. At^q already has a 
tremendous number of different isomers, and a much more complicated form of the DOZ 
arises (see Fig. 4 (b)). The DOZ cuts the complex temperature plane into three regions 
with two transition lines approaching the real axis. Comparing with the literature the 
region below 31 K can be identified as the solid phase and the region above 45 K as a fluid 
phase. Because our MC simulations are performed at zero pressure at this temperature also 
the evaporation of atoms from the cluster starts which corresponds to the onset of the gas 
phase. The phase between these two transition lines is commonly interpreted as the melting, 
isomer hopping, or coexistance region. The analysis of the order of the phase transitions is 
quite difficult in this case and will be investigated in a more systematic study. Nevertheless 
the DOZ displays in a distinct manner the phase separation for At^q and can be viewed as 
a unique fingerprint. 

In conclusion we have found that the DOZ of the canonical partition function can be 
used to classify phase transitions in finite systems. The DOZ of a specific system acts like 



a unique fingerprint. The classification scheme given above is equivalent to that given by 
Grossmann et al. but extended to the region of finite particle numbers. We have found 
that the zeros of the partition function act like boundary posts between different phases in 
the complex temperature plane. The finite size effects for the Bose-Einstein condensation 
are reflected by an dependence of the parameter ti and only a slight change of the 
parameter a which indicates the order of the phase transition. For Ar clusters the DOZ leads 
to enlightening pictures of the complex process of melting or isomer hopping, identifying in 
a distinct manner two critical temperatures supporting an old assumption of Berry et al. 
T8| , p^ . This classification scheme developed for the canonical ensemble should also hold for 



other ensembles, i.e. different experimental conditions should not influence the nature of 
the systems although e.g. the shape of the caloric curve may significantly differ. 
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FIG. 1. Schematic plot of the DOZ illustrating the definition of the classification parameters 
given in the text. 
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FIG. 2. Contour plots the ground state occupation number | rjo \ /N in the complex tempera- 
ture plane for 40, 120, and 300 particles in a 3-dimensional isotropic trap. The black spots indicate 
the locations of zeros of the partition function. 
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FIG. 3. Plots of the classification parameters a, 7, and ti versus the number of particles for a 
3-dimensional harmonic trap. 
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FIG. 4. Gontour plots of the specific heat | cy | for Arg and Arso clusters. 
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